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SUMMARY 


The method presented in NACA TN 3910 for calculating the effect of 
spanwise variations in gust intensity on the statistical characteristics 
of the response of. an airplane to continuous random atmospheric turbulence 
is used herein to calculate the effect of these variations on the lift 
directly due to turbulence. Both the horizontal and the vertical com- 
ponents of turbulence for both swept and unswept wings are considered. 
Several analytic approximations to the correlation functions and power 
spectra of atmospheric turbulence and several spanwise weighting functions 
(span loadings) are used in these calculations. 

The averaging effect of the span on the lift is shown to be very 
similar for the normal and the longitudinal components of turbulence, for 
the various span loadings and turbulence spectra considered herein, and 
for various angles of sweepback. 


INTRODUCTION 


The local fluctuations of air velocity experienced by an airplane 
flying through atmospheric turbulence constitute a random process. If 
the fluctuations are assumed to be stationary in a statistical sense, 
that is, if the statistical characteristics of the turbulence are assumed 
to be substantially invariant along the flight path, the techniques of 
generalized harmonic analysis may be used to calculate the airplane 
responses. (See ref. 1, for instance.) 

In this approach the assumption has generally been made implicitly 
that the gust intensity at any one instant is substantially the same at 
all points on the wing. When spanwise variations in gust intensity are 
to be taken into account, the statistical problem becomes more difficult 
inasmuch as the input process is now multidimensional and the conventional 
techniques of generalized harmonic analysis can no longer be used directly. 

This problem has been treated in references 2, 3* and k, and the 
results of some calculations were presented in references 2 and 3* These 
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results Indicate the magnitude of the effect of spanwise variations in 
instantaneous gust intensity on same airplane responses , including the 
lift. However , these results were very limited in scope inasmuch as 
only one atmospheric spectrum and only one spanwise weighting function 
were considered. Also, only spanwise variations in. the normal com- 
ponent of atmospheric turbulence were considered. 

The purpose of this paper is to present the results of calculations 
of the lift spectra and correlation functions for several spanwise 
weighting functions and various atmospheric spectra, both for the nor- 
mal and the longitudinal components of turbulence . To this end, the basic 
relations for the lift correlation function and spectrum derived in ref- 
erence 2 are restated and extended to the case of the longitudinal com- 
ponent of turbulence. Several analytical expressions for the spectra 
and weighting functions which have some of the desired characteristics 
are then used to perform the calculations. 

The turbulence is assumed to be stationary, homogeneous, and iso- 
tropic (although the somewhat less restriptiye assumption of axisymmetry 
would suffice for almost all the calculations presented herein, as in 
refs. 2 and 3) • Thus, the statistical characteristics of the turbulence 
is assumed to be invariant under a displacement of the time or space 
origins, as well as under a rotation or inversion of the coordinate axes. 


SYMBOLS 


b 

C 

°Ia 

c l 

c 

c 

E 

h 

^1 


span 

arbitrary parameter in the point X one-dimensional) correlation 
functions and power spectra bf case 6 

wing lift -curve slope 
section lift coefficient 
local chord 

average chord, s/b .. 

complete elliptic integral of the second kind 

dimensional unsteady-lift influence function for indicial 
gust penetration 

dimensional unsteady -lift function for ind icia] gust 
penetration 
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KLq 

k 2 

k' 

L 

L* 

L+ 

q. 

s 

r 

t 

u 

u 

v 

w 

x 

y 

y* 


d im ensional unsteady-lift function for sinusoidal gust 
penetration 

Bessel functions of the first kind, order 0 and 1 

complete elliptic integral of the first kind 

modified Bessel functions of the second kind, order 0 and 1 

incomplete modified Bessel functions of the second kind, 
order 0 and 1 

integral of Kq 

dimensionless unsteady-lift function for sharp-edged gust 
penetration 

dimensionless frequency of wave number, aiL*/u 
lift 

longitudinal integral scale of turbulence 
lateral integral scale of turbulence 
dynamic pressure 
wing area 

dimensionless lateral displacement, t]/l* 
time 

forward speed of airplane 

longitudinal (horizontal) component of turbulence 
lateral (side) component of turbulence 
normal (vertical) component of turbulence 
coordinate along the mean flight path 
lateral coordinate 


dimensionless lateral coordinate, 

b/2 



NA.CA TN 5920 


angle of attack, radians 
span ratio, b/L* 

spanwise lift influence function (lift distribution for unit 
angle of attack in reverse flow) 

autoconvolution of 7 

Fourier transform of 7 

lateral displacement 

T1 

dimensionless lateral displacement , — •— 

b/2 

longitudinal displacement 

dimensionless longitudinal displacement, £/L* 
angle of sveepback, deg 

time displacement, argument of time correlation function 
point (one -dimensional) power spectrum 

single Fourier transform of tWo-dimensional correlation 

function (for isotropic turbulence, unless stated otherwise) 

two-dimensional power spectrum, double Fourier transform of 
two-dimensional correlation function 

two -dimensional, power spectrum, double Fourier transform of _ 
two-dimensional correlation ifunction (for isotropic 
turbulence) 

dimensionless unsteady-lift function for sinusoidal gusts 
(Sears function) i 

point (one-dimensional) correlation function 

two-dimensional correlation function 

frequency of oscillation ; 

oib/2U 


dimensionless frequency. 
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Subscripts : 

e effective (averaged) 

L lift 

u, v, w components of gust velocity 


BASIC RELATIONS FOR THE CORRELATION FUNCTIONS 
AND SPECTRA OF THE LIFT 


Normal Component of Turbulence 

The material presented in this section is substantially identical 
to the treatment of this subject in references 2 and 3 and is given here 
primarily as a basis for the analysis of the longitudinal component of 
turbulence in the following section. 

Unswept wings .- The instantaneous value of the lift due to the nor- 
mal component of turbulence can be written in terms of an indicial- 
response influence function h(t,y) as 

P<» pb/2 

L(t) = /. / 

U -m b/2 

As pointed out in reference 2 , if the wing is unswept, the function h(t,y) 
can be written as the product of a function which depends only on time and 
a function which depends only on the distance along the span; that is. 


h(t,y) v(u(t-ti) ,y) dy dt x 


(1) 


h(t,y) = i h^t) y(y) 


( 2 ) 


where 


h x (t) 


C ^9S d 
U dt 




r(y) = 


Mi 
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and k 2 is the unsteady-lift function due to gust penetration normalized 

with the steady-state value. The superscript R refers to reverse flow 
and can "be disregarded for unswept wings, ‘because the lift distribution 
of an unswept wing is the same in direct and reverse flow within the 
limitations of linearized lifting- surface ‘theory. 

Equation (l) can therefore be writteh as 


L(t) = J h^tj dt-L i J ^ r(y) w(u(t-tj ,yj dy 

Hence, the correlation function for the lift can be written as 
V T ) = f_" 11 Oh) h ( t 2) dt i dt 2 ^ f ,/ /2 r(yi) r(ya) ¥ w (u(r+ t 1 -t e ),y 2 -y 1 ) ay x dy 2 ( 5 ) 

where \jr w is a two-dimensional correlation function defined by 


= w(x,y) w(ix+|,y+ri) 

where the bar designates an average, which may be performed over time or 
space, in view of the as suned stationarity and homogeneity. For turbu- 
lence which is axisymmetric or isotropic; this double correlation func- 
tion can be expressed in terms of the point (one -dimensional) correlation 
function as 


t w ( s,n) = t w (j|i^i 2 ) (M 

If the effect of the spanwise variation in gust intensity is not 
taken into account, the double integral over the span in equation (5) 

is replaced by b 2 (u (r+ti-t 2 )) • Consequently, this integral repre- 
sents, essentially, an "averaged” correlation function for the normal 
component of velocity which may be' used in the place of in 

order to take into account spanwise variations in gust intensity. 

I 

With the use of equation (4) this function can be defined as 
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and hence 


= f f h ( t l) h ( t 2) ♦w e ( u ( T+t l" t 2)) dt l dt 2 


( 6 ) 


The problem of obtaining from i|r w is identical to the one for 

the case where spanwise variations in gust intensity are ignored and, 
hence, will not be considered herein; instead, attention will be confined 
to the function ijr w (£). 


The double integral of equation ( 5 ) can be reduced to a single inte- 
gral by introducing the function T(tj) defined as the autoconvolution 
of 7 ( 7)3 namely 

p (b/2)-n 

r(n) = f / 7(7) r(y+ti) ay 

D j -b/2 

In terms of this function, equation (5) can be written as 

*We (g) = i / Q b r (n) t w (fi w) dri (7) 

In view of the symmetry of 7(7)3 the function T(ti) can also be 
written as 


, p( b/2)-ii 

nn). = £ / 7(7) 7 (y+n) &y 

b J - tj/2 


and it must satisfy the property 


( 8 ) 



di| = b 


(9) 


* The power spectrum of the lift can be obtained by calculating the 
Fourier transform of the correlation function given by equation ( 3 ); 
namely. 
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<pi>) = i r e "‘ icurr t L (T) dT 

O —oo ; 

= |h-l(<d )| 2 cp Ve (qi) (io) 

where H^cd) is the Fourier transform of h^(t) and cp We (aj) is an 

"averaged" spectrum which bears the same relation to the point spectrum 
as the correlation function does to : and which is defined by 

'PwpM r e ’ 1CDT t v „(U T ) dt 

e * J -oo e 

p oo -i^| 

- £ L e V‘> a6 (11) 

The function H]_(ui) represents the ponrplex amplitude of the lift 

response to sinusoidal gusts which are uniform along the Bpan and can 
be expressed as 


HiM = ^ 0(§) 

where 0 represents the unsteady-lift function for sinusoidal gusts, 
sometimes referred to as the Sears function. ^ As in the case of the 
correlation function, attention will be cbnfined to q> w because the 

manner of obtaining the lift spectrum from it is straightforward. (See 
eq. (10).) , 

As pointed out in references 2 and 3; an alternative method of 
obtaining cp We consists in using two-dimensional spectra cp w (m>,T]) and 

cp w (o>]_,G^ defined for a general two-dime'nsional correlation function 

as ! ■' 


?w(^Tl) 



? w (5,ti) d| 


(12) 



MCA Off 5920 


9 


and 


s ^ 


/: r 




d s d y 


.00 -CO 


iaio 

i r 00 ~r T 


(15) 


f e U $vr(“l^) d y 


For axisymmetric or isotropic turbulence. 


?V(^Tl) - ^ y' e ^ t w (fe^n 2 ) U 


and ^ w (o>L,a>2) becomes a function cp^. of only one variable, 

cd s \jay]_^ + o>>2; this function can be obtained from i|r w directly as 
follows : 

- J? f 0 iJ o($ t) t„(5) « 

The spectrum cp^Ccs) can then be obtained directly from the fol- 
lowing expressions, instead of using equation (ll) : 

r* "k 

cpw e (co) = i J r(rj ) c^Ccjd,!}) dll (HO 

or 

9v e ^“) “ J ? 2 K) Tw^a^-Ki^ 2 ) dn^ (15) 

where r(^) is the autoconvolution function defined in equation (8), and 

A 

r(cn) -is the Fourier transform of 7(y), namely 
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- * c 


( 16 ) 


This function is real and symmetric "because ?(y) is symmetric. 
Swept wings . - For a swept wing, equation (l) is replaced by 

L(t) = J' S / li ( t i> y ) w ( u ( t "^i) + |y| ta ^#y) <2y cl^i 


( 17 ) 


Whether the function h(t,y) can still be written in the form of equa- 
tion (2) is not known at present. However, in the following equation the 
assumption is made that equation (2) is still valid. The averaged corre- 
lation function i!r w (£) can then be written'as 

e I 

y ( y :0 y ( y s) tv (^ + (| y 2| -Jyij) 2 + (y 2 -yi) 2 j ^1 ^2 ds> 

Inasmuch as y 2 and y-j_ no longer Appear only as the difference 
(y 2 - y-^ , equation (l 8 ) cannot be simplified to the form indicated by 
equation ( 7 ) and must in general be evaluated numerically. 

The spectrum cp Wg (m) can then be ob-jtained by calculating (numerically) 
the Fourier transform of ® ie 3i rec i expressions for cp Wg given in 

equations (l4) and ( 15 ) are not valid for: the swept wing. 


Longitudinal Component of Turbulence 

Unswept wings .- The instantaneous vaJLue of the lift due to the longi- 
tudinal component of turbulence can be written as 


L(t) = 2a J' J h(ti,y) u^J (t-tq) ,y) dy dt] 


(19) 


where the function h(t,y) is the same iqdicial-response influence func- 
tion as the one used in equation (l) , provided that the angle of attack a 
is constant along the span. 

The correlation function for this lift can be written in the form 
of equation ( 3 ) • However, equation (4) is not valid for the longitudinal 
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component . Instead, for axisymmetrical turbulence the two-dimensional 
correlation function for the longitudinal component can he written in 
terms of the point (one -dimensional) correlation function for the lon- 
gitudinal and lateral components as follows : 

*u<^) = t u (^ 2 +tl 2 ) + -~—z t v (l|sV) (20) 


|2 + tj2 


i + ti c 


If the function on the right side of equation ( 20 ) is substituted 

for in equations (5) and (j) } corresponding expressions for 

the averaged correlation function iJr We (|) are obtained in terms of the 
same functions y(j) and r(q) used previously, namely 


W S) ' S Jo ™ PT^ ^(1^) + 


1 Mr 

„p p 

%r + tt 




dri 


Similarly, if the same substitution is effected in equation (12), an 
expression for cp u (cn,q) is obtained; namely. 




+ T] 


2 ^v 


■(^ 2 +n 2 ) 


d| 


With this expression and the function r(t]) used previously, an averaged 
spectrum 7^(01) can then be obtained from equation (ll) by use of the 

function V (|) or, also, from equation (l 4 ) by use of 9 u (a),T]); thus. 


c Pu e ( m ) = ^ / r(n) 9 u (co,Ti) dT] 


Although a function <p u ^a>j_,cD2) can be calculated analogous to the func- 
tion cp w (o>]_,ui2^, this function cannot be expressed in terms of ^o>j_2 4. ^2 
alone; therefore the advantages of this approach are largely lost. 

From this spectrum (a>) , the lift due to the longitudinal component 
of turbulence can then be obtained by use of the relation 
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CPlCoj) = lkt 2 !%(«>) | 2 cp^Coi) 


( 21 ) 


As pointed out in reference 2, the lifts due to the normal and longi- 
tudinal components of turbulence are statistically independent for iso- 
tropic turbulence; consequently, the spectrum of the total lift due to 
both components is the sum of the spectra given by equations (lO) and (21) 


Swept vings .- The analysis of the lift on a swept wing due to the 
normal component of turbulence can also be modified to apply to the 
longitudinal component in a straightforward fashion . The function Tjr w 

under the integral in equation (l8) must : be replaced by 


C e 1 (N ~ Ini ) tan A ] 2 K(s) , >2 - n) 


R c 


* V (R) 


where 


R = ^ i + ^|y 2 | - |yi|)tan Ij + (y 2 - yi) : 


and the evaluation of the resulting double integral for t\i e (l) must 
again be effected by numerical methods, As must the calculation of 3*^(0}) 
from this function W&). 


If, in the preceding section (dealing with the normal component 
of turbulence for swept and unswept wings) the turbulence is assumed to 
be only axisymmetric and not isotropic, the functions tvU) and *w(fc> 
are not identical; the relation between them depends on the variation of 
the statistical characteristics of the turbulence in a direction perpen- 
dicular to the xy -plane. However, if coipplete isotropy is assumed, T|r v 
and \|r w are identical, and a definite relation exists between t1t v , i|r w , 
and which will be given in the following section; equation (10) can 

then be expressed in terms of \|r u alone! or f w alone. However, little 
practical advantage is gained in this manner. 


WEIGHTING FUNCTIONS AND SPECTRA USED. IN THE CALCULATIONS 


Weighting Functions 

The lift influence function y(y) , which has been Identified with 
the lift distribution due to ■unit angle Of attack in reverse flow, serves 
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as a weighting function in the integral expression for T(r We , 

and cp^ in terms of the point correlation functions and spectra, The 

functions used herein for the purpose of indicating the effect of changes 
in y(y) on the averaged coir elation functions and spectra are given by 
the simple analytic expressions listed in table 1 and are shown in fig- 
ure 1; they will' be referred to as rectangular, elliptic, triangular, and 
parabolic distribution, respectively. All functions satisfy the condition 

f b/2 

/ ny) dy = b 
J -b/2 

which follows from the definition of y(y) . The function designated as 
rectangular distribution is the one that has been used in references 2 
and 3- The corresponding functions r(i}) and P(cd) are also given in 
table 1, and the function r(ri) is shown in figure 1. . 

For functions other than those considered herein (which can be 
obtained for any given wing from reference 5 or a similar compilation 
of spanwise lift distributions), the functions r(rj) and f ( cd) must be 
calculated numerically. The integrating factors given in table 2 

have been found useful for calculating T(ti). These factors apply to 
the case where values of y(y) are used, at every tenth of the semispan 
and where r(T|) is obtained at every tenth of the span by use of the 
relation 


.. 9-2i 

r (ii) - b 11 ’'(V'i) 

with 

T). = ib yj = -sL — 

1 10 ^ 10 2 

where i = 0, 1 . . .9 and j = -9, -8 . . .8 ,9* The value of r(b) 
is always zero; therefore no factors are listed for i = 10. These fac- 
tors apply to functions /(y ) which go to zero at the wing tips as the 
square root of the distance from the tips, as spanwise lift distributions 
in subsonic and supersonic flow usually do. 


Point Correlation Functions and Spectra 

General relations .- In isotropic turbulence the mean-square values 
of the u, v, and w components of the turbulence are identical; that 
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is, = v2 = The point correlation' functions for these components 

(for flight parallel to the u component) can "be expressed in terms of 
each other by means of the relation given; in reference 6; namely. 


^ v (!) = t w (!) 


+v(D = V U (S> + |^ 1fu(l) 


( 22 ) . 


This ordinary differential equation for i^ u (|) with the boundary condi- 
tion t u (o) 4 t w (o) can be solved to yield the following expression: 


* u (s) - 



(25) 


Hence, the power spectra corresponding to these correlation func- 
tions must satisfy the relation 



(24) 


Moreover, the solution of this equation for (^(o) with the boundary 
condition cp u (°°) = 0 yields 


<P u (o>) 


= 2cu 



<pw,H) 




dto^ 


(25) ’ 


The Integral scales of turbulence L* 


and are defined as 


L* =s 
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and 


r“ 

L + =~~ / t w (l) U 

w 2 J 0 

- |^(o) 

v 2 


As a result of equation (22) 



for isotropic turbulence, as can be shown by integrating the second term 
on the right side of equation (22) by parts. Also, aB a result of equa- 
tion (22) , the moment of the lateral correlation function can be shown 
to be zero; that is, 



I d| - 0 


provided that ^ w (£ ) goes to zero more rapidly than l/g 2 when | 

becomes very large - a condition which any actual correlation function 
must satisfy. This condition on the moment of Tjr w (|) implies that for 
isotropic turbulence any actual lateral correlation function must be 
negative for some values of 5 . 

Specific expressions used in calculations .- The point correlation 
functions and spectra used in the calculations are listed in table 3* 

The one labeled "case l" is the one used previously in references 2 
and 3 • Although the spectra do not have finite moments (so that the 
mean-square values of the derivatives of u, v, and w are infinite) , 
they appear to fit available experimental data on turbulence in the atmos- 
phere and in wind tunnels fairly well, except at very high and very low 
frequencies. At very low frequencies some of the responses of airplanes 
are very small; thus, despite their theoretical shortcomings, these 
spectra are often very useful. 

The functions of case 2 are based on an exponential approximation 
to the lateral correlation function and suffer not only from the short- 
comings of the functions of case 1 but also violate the condition of zero 
moment for the function i|r v (| ) . However, case 2 represents some results 
of turbulence measurements in wind tunnels fairly well. (See ref. J, 
for instance.) 
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The functions designated by case 3 and case 4 represent fittings of 




so that 


t u (l) and Tlr w (t), respectively, by functions of the form 

-K'm^ 

the spectra contain the factor e . Although the functions of case in- 

violate the condition of zero moment for : ^r w ( £ ) , both cases may be used 
to represent all or part of some wind-tunnel measurements of turbulence. 
(See ref. 8.) Furthermore, turbulence theory (ref. 9, for instance) indi- 
cates that a certain frequency range of the turbulence spectrum should con- 
tain a factor e" K a)2 . , 

The spectrum cp u (co) of case 5 represents, essentially, an interpola- 
tion formula suggested in reference 10, where for a certain frequency 
range, the three-dimensional (Heisenberg) spectrum of isotropic turbulence 

r - 2 ”! “IT /6 

is assumed to be proportional to coMl + ^ so that the corre- 

sponding one -dimensional spectrum cp u (cb) ; used herein would be proportional 


to 


1 + 


fV 


-5/6 


The constant of proportionality can be determined 


*u<°) - 


u 2 , 


from the fact that 
the scale of turbulence, 
one given in table 3 ; the other expressions then follow from equation (22) 
and the definitions of the spectra. The ‘constant 0.59253 represents 

and the constant 1-339 represents jr(^] } where r repre- 


and the constant cjcq . can be related to 
The resulting expression for ilr u (|) is the 




r V 

J)’ — — 

sents the gamma function. These correlation functions violate the condi- 
tions on the derivatives of a correlation function. or the moments of its 
spectrum, just as the ones of case 1 do; but for a certain range of dis- 
tances or frequencies, they may serve as. useful approximations in some 
instances . 


Case 6 represents case 1 with a correction which makes the second 
moment of the spectra finite; consequently, the mean-square values of 
the first derivatives exist, although thp higher moments, and, hence, the 
mean-square values of the higher derivatives' do not exist. The con- 
stant C is arbitrary and can be 1 used tb select a frequency at which 
the point spectra start deviating from a: second -power variation with 
frequency and decreasing more rapidly, bn the basis of the measurements 
reported in reference 1, a value of Q =• 50 appears appropriate and was 
used in the calculations of the present paper . The expressions given in 
table 3 for case .6 are not entirely consistent with those for the other 
cases. In order to be consistent, a should be multiplied by (l + C”^) , 
k 1 should be divided by this quantity, and the expressions given for 
cp u (a>) and cp w (m) should be divided by this quantity. However, by virtue 
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of the fact that this quantity differs from unity by only O.Ol percent 
for C = 50, this refinement has been ignored in order to simplify the 
expressions ; it might have been worthwhile if a smaller value of G had 
been selected. 


The normal point correlation functions considered herein are shown 
in figure 2, the longitudinal point correlation functions in figure 3, 
the no rmal point spectra in figure b } and the longitudinal spectra in 
figure 5 . In a] 1 these figures, the point correlation functions and 
spectra are given by the curves for (3 = 0 when not specifically referred 
to as point correlation functions and spectra. The correlation functions 
for case 6 with C = 50 are virtually indistinguishable from those for 
case 1, except for very small values of | ; that is, |/L* less than 0.1. 
Similarly, the spectra are indistinguishable except for very high fre- 
quencies (cuL*/u greater than 50 ), for which the spectra of case 1 
attenuate as •ar*2 whereas those of case 6 attenuate as crr^. 


Which of these cases best represents the actual characteristics of 
atmospheric turbulence is not known at present. Each case (except pos- 
sibly cases 2 and *0 has some features which atmospheric turbulence is 
likely to possess, but each also has certain shortcomings. Case 1, by 
virtue of its simplicity and the fact that measured atmospheric spectra 
have large ranges in which they appear to attenuate as or 2 , is probably 
the best choice for calculations where the behavior of the spectra at 
very high frequencies is not very important; when it is important, the 
correction represented by case 6 is likely to be sufficient for many 
purposes to take the more rapid rate of decrease of the spectra at high 
frequencies into account. 


RESULTS AND DISCUSSION 


Two-Dimensional Spectra . . ■ 

The two-dimensional spectra 9 ^( 03 , 11 ), 9 ^( 03 , 11 ) , and cpy-Ccu) are 

given in table 4 for five of the six point correlation functions and 
spectra given in table 3* Again, for case 6 the values of k’ .. and of 
the spectra (^((JOjT]) and 9 ^( 03 , T]) should be divided by 1 +. C“ 2 , and 

the value of 9 ^ 03 ) should be divided by the square of this quantity 
if this refinement is considered worthwhile. 
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Mean -Square Value of the Averaged Normal and Longitudinal 
Component of Turbulence 

The mean -square values of the averaged normal and longitudinal com- 
ponent of turbulence defined by 

w e 2 = ^w e (°) 

and 

Ug 2 a tug(°) 

have no direct physical significance. If the wing responds equally to 
sinusoidal gusts of all frequencies , so that 0 (<x>) = 1 , the mean-square 
lift is proportional to these quantities ; although for actual unsteady- 
lift functions 0(o>) they have no direct relation to the lift. Nonethe- 
less , they serve as convenient quantities to normalize the effective or 
averaged correlation functions and will be used for that purpose. The 
power spectra can also be normalized with them, but the point mean-square 
value serves equally well for this purpose.. 

For isotropic turbulence the two quantities w e 2 and Ue 2 can 
readily he shown to be identical for unswept wings; therefore a distinc- 
tion is made between these quantities in! this section for swept wings 
only. 


The values of w e 2 calculated for the various point spectra are 

shown in figure 6 (a) for unswept wings and rectangular loading as a 
function of the span ratio b/L*. The analytical expressions for these 

curves are given in table 5 ■ The manner. in which w e 2 decreases with 
increasing span may be seen to depend to 'a large extent on the nature of 

the point correlation functions for small span ratios (where w e 2 differs 

relatively little from w 2 ) . For very large spans (that is, as 0 -*«>), 
the following asymptotic relation can be shown to exist: 



(26) 


Therefore the curves shown in figure 6 (b), with the exception of the one 
labeled "elliptic loading," would tend to coalesce for very large span 
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ratios. This region of span ratios is not of direct interest for air- 
planes, but may be of interest for tests of the response of airplane 
models to wind-tunnel turbulence. 

For a swept wing, Ug^ and w e ^ are not Identical. They have been 
calculated (fig. 6(a)) for rectangular loading and the point correlation 
function of ease 1 for a span ratio b/L*cos A = 0.5 as a function of A. 
(The angle A is there designated as the angle of sweepbaek, but could be 
the angle of sweepf orward as well . ) They may be seen to differ relatively 
little from each other or from their value at zero sweepbaek. 

The effect of the shape of the loading function y(y) on w e 2 may 

be seen in figure 6 (d), where Wg^/w^ is shown as a function of b/L* 
for unswept wings, for the point correlation function of case 1 , and for 
four different loading functions . For small span ratios the effect of 
changes in 7 is fairly small; only for very large span ratios does the 
effect become significant, as may also be seen from equation (2 6 ) . Simi- 
lar calculations have also been made numerically for the point correla- 
tion function of case 5 for rectangular and elliptic loadings; the results 
are shown in figure 6(b) and corroborate the conclusion reached from the 
results obtained with the point correlation functions of ease 1. 


Correlation Functions of the Averaged Components 

of Turbulence 

The functions i|r We and calculated for the point correlations 

of cases 1 to 5 for unswept wings with rectangular loadings are shown in 

figures 2 and 3 , where they have been normalized with v e ^. The corre- 
sponding analytical expressions are given in tables 5 and 6. Inasmuch as 
no tables of incomplete modified Bessel functions were available, the func- 
tions for cases 1 and 5 were calculated numerically and those for case 2 
were obtained from reference 7 where they were also calculated numerically. 

The effect of the span ratio b/L* on the normalized correlation 
functio ns may be seen to be fairly small and in some cases nil, namely 
for case 3 and i|r w for case 4. This effect appears to be sub- 

stanti ally the same on both tug and - as “ay be seen by comparing 

figures 2(a) and 3(b). The nature of the 'effect depends on the point 
correlation functions. For instance, in the range of |/L* from 0 to 2, 
the normalized functions i(fw e tend to increase with p for cases 1, 2, 

5, and 6, but to decrease or be unchanged for cases 3 and h. The func- 
tions labeled P-> » pertain to span ratios so large that w e ^ is given 
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by equation (26) with sufficient accuracy and are, again, not directly 
pertinent to actual airplanes but are shown primarily as a theoretical 
limiting case . 

The effective correlation functions and for swept wings 

with b/L*eos A = 0.5 and rectangular loading and for the point corre- 
lation functions of case 1 are shown in figures 7(a) and 7(^)- When 

these functions are normalized with w e 2 and u^ 2 , respectively, they 
are not affected by sweep to an appreciable extent. These effective 
correlative functions were calculated from the integral expressions 
given in this paper by use of a combination of numerical and graphical 
techniques . 

Some calculations of. the effective correlation functions were also 
made for elliptic span loading or weighting functions y( y), by use of 
numerical methods and are shown in figures 2(a), 2(d), and 3(b). The 
differences between them and those for rectangular loading may be seen 

to be very small if the functions are normalized with w e ^. 


Spectra of the Averaged Components of Turbulence 

The functions cp Wg and cp^ calculated for the various point corre- 
lation functions for unswept wings with Rectangular loadings are shown in 
figures 4 and 5* In figures 4(b), 4(g), and 5(c), these functions have 
been normalized with w 2 , b ut i n the other parts of these figures they 
have been normalized with w e 2 and thus ; represent directly t he F ourier 
transforms of the functions tw e and normalized with w e 2. The 

analytic expressions for these functions; are given in tables 5 and 6. 

An examination of these figures indicates that, regardless of whether 

• — « - 

the spectra are normalized with w e 2 or W 2 , the effect of spanwlse vari- 
ations in gust intensity on the spectra appears to be larger than the 
effect on the correlation functions. From figures 4(b) and 5( a ) may be 
seen the fact that this averaging process tends to reduce the magnitude 
of the spectra at all frequencies, but to the largest extent at high 
frequencies, as may be expected. For the point spectra of cases 1 , 5> 

and 6, which attenuate at high frequencies as co -2 , cd“5/3 } and cd"^ 

respectively, the effective spectra attenuate as cu~3 } oj- 8/5, and txT^ , 

respectively. This fact may be seen from figures 4 and 5 and follows 
from the asymptotic relation for 

<Pw e (c») ~ f 


(27) 
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and from a similar relation for for instance, for case 1 and rec- 

tangular loading 


<Pv e (“) 3 

v2l*/u 3 k ’ 5 

whereas 

qvk°) 3 

- - - 

w 2l*/u ' 

Again, the nature of the effect of averaging is substantially the 
same for cp^ and qp^, but appears to depend to some extent on the 

point correlation function, the effect being fairly large for the point 
correlation functions of cases 1 , 2 , 5 * and 6 and relatively small or 
nil for those of cases 3 and. 4. 

The effect of sweep on the effective normal and longitudinal spectra 
may be seen from figures 7 (c) and 7(3.) to be very small for the given 

value b/L*cos A = 0.5* provided that the spectra are normalized with w e 2 

and Uq 2, respectively. These spectra were calculated by performing the 

Fourier transforms of and V for these cases numerically, and the 

results of these calculations are somewhat uncertain at high frequencies. 
Therefore these functions are shown only for frequencies up to cuL*/U = 20 
and may be valid only for oiL*/U below about 10. The possibility exists, 
therefore, that the spectra for swept wings differ from th ose f or unswept 

wings at very high frequencies, even when normalized with w g 2 . 

The effect of changes in the shape of y(y) on the effective spectra 
may be seen from figure 4(g), where the spectra are shown for unswept wings 
with b/L* = 0.5 and for the point correlation functions of case 1 , but 
for several loadings . These spectra were normalized with w^ in order to 
emphasize the effects of the’ shape of the loading. Even so, these effects 
are small, except at very high frequencies, as may als o be seen from equa- 
tion ( 27 ) • If the functions had been normalized with w e 2 , these effects 
would have been even less pronounced. For some other point correlation 
functions, calculations have been made numerically for ellip tic as well 
as rectangular loadings. The results, normalized with w g 2, are shown 
in figures 4(a), 4(e), and 5 (b) and may be seen to differ very little 
from the corresponding ones for rectangular loading. 



22 


NACA TN 3920 


Effective Scales of Turbulence 

Effective or apparent scales of turbulence can be defined in a man- 
ner analogous to that for the actual scales of turbulence, namely 


and 



These parameters may therefore be deduced from the relations for the 
spectra given in tables 5 and 6 for <jo = 0 and are shown in figure 8 
as functions of the span ratio b/L* for "the various point correlation 
functions used herein. Figures 8(b) and 8(d) jaertain to unswept wings, 
and figures 8(a) and 8(c) pertain to swept [wings with b/L*cos A = 0.5 
and rectangular loading. I 

The averaging effect of the span on these_ scales depends on the 
point correlation functions and appears to be fairly small for small 
span ratios; it is similar for both L* e hnd L^ e * For instance, L^ e 

tends to increase as the span increases for the point correlation func- 
tions of cases 1, 2, and 6 and to decrease for case 3# it is inde- 
pendent of the span, for case 1. 

As may be seen from figures 8(a) and 8(c) the effect of sweep on 
L^ e and L* e appears to be very small fob given values of b/L*cos A. 
Similarly, the effect of changes in the shhpe of y(y) on L^ e and L* e 
may be seen to be fairly small from figured 8(b) and 8(d), where for some 
point correlation functions results are shower hot only for rectangular 
loadings but also for elliptic, parabolic,! and^ triangular loadings. 
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In figures 8(1 d) and 8(d) the rays for constant hj L"^ e and b^L* e 
have been drawn in. The purpose of this addition to the figures was to 
facilitate the problem of determining a value of L* or L"^ from lift 
measurements on a wing of finite span. In this problem cp We and <p Ue 

can be deduced if the unsteady-lift function and the lift spectrum are 
known. Hence, L* e or L^ e is known, as is the span of the wing; 
therefore, the ratio hjl,* e or b^L^„ is known. The intersection of 

the ray for that ratio with the curve for the case for which the shapes 
of the averaged spectra most nearly resemble that of the spectrum obtained 

from the measurements establishes the values of or L* e jh* and 

of b/L*, from any of which the value of L* can be determined directly. 

CONCLUSIONS 


The correlation functions and spectra for a suitably weighted 
instantaneous average of the intensity of the normal and longitudinal 
components of atmospheric turbulence have been calculated for isotropic 
turbulence by the method of MCA TN 3910. Several analytic expressions 
were used for the point correlation functions or spectra and for the 
lift influence function (weighting function) . These averaged correla- 
tion functions and spectra bear the same relation to those for the lift 
as do the point correlation functions and spectra, except that they take 
into account the instantaneous variations in gust intensity over the 
span; whereas the use of the point correlation functions and spectra 
implies that the gust intensity is assumed to be constant along the span 
at any Instant. Consequently, the following conclusions concerning the 
averaged or effective gust intensity pertain directly to the power spec- 
trum of the lift: 

1. The effect of spanwise variations in gust intensity is to decrease, 
essentially, the mean-square gust intensity experienced by the wing. The 
extent of this decrease is small when the span is small, but increases as 
the span increases ; the rate at which it increases depends on the shapes 
of the point correlation functions and of the weighting functions. The 
effective mean-square longitudinal and normal components are the same for 
unswept wings (in isotropic turbulence) but differ slightly for swept 
wings. 

2. Apart from this overall decrease in level, the correlation func- 
tions are relatively insensitive to changes in span, provided that the 
span is small compared with the scale of turbulence as is the case for 
actual airplanes. The small changes in the shape of the correlation 
functions depend primarily on the shape of the point correlation function 
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and, to a lesser extent, on the shape of! the weighting functions and the 
angle of-sweepback, and are similar for the normal and longitudinal com- 
ponent 'of turbulence. 

3* The instantaneous spanwise variations in gust intensity result 
in a decrease of the spectra at all frequencies, but most significantly 
at the highest frequencies . The extent of this decrease becomes greater 
as the span increases and depends on the ! shape of the point power spectra, 
but is less sensitive to variations in trie spanwise weighting function 
and to the angle of sweepback, except at 'very high frequencies; also this 
decrease is similar for the normal and longitudinal component of turbulence . 

4. The scales of turbulence are affected relatively little by changes 
in span, at least for spans which are smhll compared with the scale of 
turbulence . The slight changes in the scales depend primarily on the 
point correlation functions or spectra and may represent either an increase 
or a decrease of the apparent scales; they are not affected appreciably by 
the shape of the spanwise loading function and the angle of sweepback. 


Langley Aeronautical Laboratory, 

National Advisory Coimnittee for Aeronautics, 
Langley Field, Va., November 2, ! 1956. 
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TABLE 1.- WEIGHTING FUNCTIONS /(y) AND CORRESPONDING FUNCTIONS r(ri) AND ?(©) 


Shape of weighting 
function 

r(y) 

r(tj) 

F(o>) 

. 

Rectangular 

1 

2 - n* 

Bin as* 
<# 

Elliptic 

& - y * 2 

+ ,*) [(1 ♦ I n* 2 ) e(^,) - o* k(|j| 

J-i(co*) 
2 -“5 ~ 

CO* 

Triangular 

2(1 - |y*j) 

< 

|(V - 6tj* 2 + 3n^) 0 ^ T] ^ | 

|(8 - i2n* + 6n*^ - n* 5 ) | ^ n < h 

n 1 - COB CO* 

* *2 
fly*** 

Parabolic 

■ 

^(32 - 40t|* 2 + 20T]* 3 - i!* 5 ) 

sin co*. - 0# cos cd* 

3 *3 
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TABLE 2.- THE INTEGRATING- FACTOR I-y FOR NUMERICAL EVALUATION OF FROM 


J 

1 = 0 

1 = 1 

1 = 2 

rc\ 

'1 

-H 

1 = 

■ 

1 = 5 

1 s 6 

i = 7 

1 = 8 

1=9 

-9 










O.15708 

-8 









0.09562 


-7 








0.13333 

■25333 


-6 







0.13333 

.12000 



-5 






0.13333 

.13333 

.30170 








0.13333 

.13333 

.22896 




-3 




0.13333 

•13333 

.26667 

.25333 




-2 



0.13333 

•13333 

.26667 

.12000 





-I 


0.13333 

.13333 

.26667 

.13333 

.30170 





0 

0.13333 

.13333 

.26667 

•13333 

.22896 






1 

.13333 

.26667 

•13333 

.26667 

.25333 






2 

.26667 

.13333 

.26667 

.12000 







3 

■13333 

.26667 

.13333 

.30170 







4 

.26667 

.13333 

.22896 








5 

.13333 

.26667 

•25333 








6 

.26667 

.12000 



■ 






7 

.13333 

.30170 









8 

.25000 










9 

.20000 
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table 3.- Kiim OTHRmriair jujctime akd hfecsba 


Case 

tj t) 
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u 2 l*/u 

$v(®) 

i?L*/U 

1 
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TABUS 3*- AVHIACSE CCRHBLAKOT HHCTinSB AM) POKER 3PBCTBA FOR VHCECCAL COMPOSHH? OP ttFBUJJJd 


V>J 
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[Rpcti ii g n l« r loading) unavopt vinga] 




EAJ3LE 6 .- AVERAGED CORRELATION FONCTIONS AND POWER SPECTRA FOR 
LONGITUDINAL COMPONENT OF TURBULENCE 
[Rectangular loading; unswept wings] 










Dimensionless lateral displacement, tj 


Figure 1.- Spartwise loading functions T'(y) and corresponding fune 
tions r(-n) used in calculations . 






firngncjoriisss Icngltydino) dispiPC6ni6nt } 


(a) Case 1. 


correlation functions for unsvept 

wings with rectangular loading (unless specified otherwise). 
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Dimensionless iongitudinai displacement , 

(c) Gases 5 end 4. 


Figure 2.- Continued 







Dimensionless longitudnal displacement 


(d) Cases 5 and 6. 
Figure 2.- Concluded. 







Dimensionless longitudinal displacement. 


L* 


(t>) Averaged correlation functions for case 1 for rectangular loading 

(unless specified, otherwise) . 


Figure 3 .- Concluded. 
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(a) Case 1. 

Figure 4 .- Point and averaged normal power spectra for unsvept wings 
with rectangular loading (unless specified otherwise) . 
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Dimensionless frequency or wave number, 


(b) Case 1; normalized with w 2 . 


Figure 4.- Continued. 





Dimensionless frequency or wave number, ^j=- 


(c) Cases 1 and 2. 
Figure 4.- Continued. 
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Dimensionless frequency or iwave number, 


(d) Cases 5 and ■ 
Figure 4.- Continued. 
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Dimensionless frequency or wave number , ‘‘jj-* 


6 8 10 20 30 40 50 


(g) Case 1; various loadings. 
Figure 4.- Concluded. 
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Dimensionless frequency : or wove number, ^ 

(a) Point spectra. 

Figure 5 • - Longitudinal power spectra for unswept wings . 
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Dimensionless frequency or wave number , 


(b) Averaged, spectra for case 1; rectangular loading ( unl ess specified 

otherwise ) . 


Figure 5.- Continued. 
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(c) Averaged spectra for case 1; rectangular loading; normalized with 


Figure 5*- Concluded. 
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(c) Unswept wings with various loadings; case 1. 
Figure 6.- Concluded. 


Dimensionless longitudinal displacement , 

(a) Normal correlation functions. 

Figure 7.~ Averaged correlation fimeiions and spectra for swept wings 
with rectangular loading and cos A = 0.5; case 1. 
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Dimensionless frequency or wave number , 
(c) Normal spectra. 


Figure 7*- Continued 




Uq M 

— w e— 008 

T u 

.006 



NACA TN 3920 


55 



(a) Normal scales for swept wings with rectangular load ing ; 


777 cos A = 0.5. 


b 

L + e 



Dimensionless span, 


(h) Normal scales for unswept wings. 
Figure 8 .- Effective scales of turbulence . 
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Angle of sweepback, A, deg 


(c) Longitudinal spectra for swept, wings with rectangular loading; 

— cos A = 0-5* 

L* 


b 



(d) Longitudinal spectra for unswept wings. 
Figure 8.- Concluded. 
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